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Abstract 

The semileptonic I? c -meson decays into the heavy quarkonia J/ip{r] c ) and a pair of 
leptons are investigated in the framework of three-point sum rules of QCD and NRQCD. 
Calculations of analytical expressions for the spectral densities of QCD and NRQCD cor- 
relators with account for the Coulomb-like a s /v terms are presented. The generalized 
relations due to the spin symmetry of NRQCD for the form factors of B c —* J/ifj(rj c )li'i 
transitions with / denoting one of the leptons e, or r, are derived at the recoil momen- 
tum close to zero. This allows one to express all NRQCD form factors through a single 
universal quantity, an analogue of Isgur-Wise function at the maximal invariant mass 
of lepton pair. The gluon condensate corrections to three-point functions are calculated 
both in full QCD in the Borel transform scheme and in NRQCD in the moment scheme. 
This enlarges the parameteric stability region of sum rule method, that makes the results 
of the approach to be more reliable. Numerical estimates of widths for the transitions 
of B c — > J '/tp(r] c )li'i are presented. 
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1 Introduction 



Recently, the CDF collaboration reported on the first experimental observation of the B c me- 
son, the heavy quarkonium with the mixed heavy flavour [Qj. This meson stands among the 
families of charmonium cc and bottominum bb in what concerns its spectroscopic properties: 
two heavy quarks move nonrelativistically, since the confinement scale, determining the pres- 
ence of light degrees of freedom (sea of gluons and quarks), is suppressed with respect to the 
heavy quark masses mq as well as the Coulomb-like exchanges result in transfers about a s rriQ, 
which is again much less than the heavy quark mass. So, the nonrelativistic picture of binding 
the quarks leads to the well-known arrangement of system levels, which is very similar for the 
families mentioned above. The calculations of fee-mass spectrum were reviewed in |2|, [§• So, 
we expect 

M Bc = 6.25 ± 0.03 GeV, 

when the measurement gave 

Mgf = 6.40 ± 0.39 ± 0.13 GeV. 

There is an essential difference in the production mechanisms of heavy quarkonia cc, bb 
and B c . To bind the b and c quarks, one has to produce four heavy quarks in the flavour 
conserving interactions^, which allows one to use the perturbative QCD, since the virtualities 
are determined by the scale of heavy quarks mass. Thus, we see that the production of B c is 
relatively suppressed a(B c ) / a(bb) ~ 10 -3 because of the additional heavy quark pair in final 
states. The basic peculiarities of production mechanisms appear due to the higher orders of 
QCD even in the leading approximation: the fragmentation regime at high transverse momenta 
much greater than the quark masses and a strong role of nonfragmentational contributions at 
p T ~ rriQ, which can be exactly calculated perturbatively; a negligible contribution of octet 
mechanism [|J] because there is no enhancement due to a lower order in a s . The predictions 
for the cross sections and distributions of B c in various interactions are discussed in []5| , where 
we see a good agreement with the meausrements of CDF [1]. 

In contrast to cc and bb states decaying due to the annihilation into the light quarks and 
gluons, the B c meson is a long lived particle, since it decays due to the weak interaction. 
The lifetime and various modes of decays were analyzed in the framework of a) potential 
models ||, b) technique of Operator Product Expansion in the effective theory of NRQCD 
, considering the series in both a small relative velocity v of heavy quarks inside the meson 
and the inverse heavy quark mass ||, c) QCD sum rules |J |TD[ applied to the three-point 
correlators [11], 0, |i~3f . The results of potential models and NRQCD are in agreement with 



each other. So, we expect that the total lifetime is equal to 

t Bc = 0.55 ± 0.15 ps, 
which agrees with the experimental value given by CDF [1]: 

t% = 0.4613:8 ± 0.03 ps, 

within the accuracy available. 

1 We do not consider for the production in weak interactions. 
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Further, the consideration of exclusive B c decays in the framework of QCD sum rules 
indicated that the role of Coulomb corrections to the bare quark loop results could be very 
important to reach the agreement with the other approaches mentioned (T^|. This requires 
working out the a s /v corrections in NRQCD, which possesses a spin symmetry providing some 
relations between the exclusive form factors. For the semileptonic decays, such a relation was 



derived in [T4[. Note, that the CDF Collaboration observed 20 events of B~£ —> ipe + (fi + )v, 
so that the consistent calculation of semileptonic decay modes is of interest. For theoretical 
reviews on the B c meson physics, see Jl5] . 

In this paper we perform a detailed analysis of semileptonic B c decays in the framework 
of sum rules in QCD and NRQCD. We recalculate the double spectral densities available 
previously in [|TT[] in full QCD for the massless leptons and add the analytical expressions for 
the form factors necessary in evaluation of decays to massive leptons and P-wave levels of 
quarkonium with different quark masses. We analyze the NRQCD sum rules for the three- 
point correlators for the first time. We derive generalized relations between the NRQCD 



form factors, which extends the consideration in |T3j , because we explore a soft limit of recoil 
momentum close to zero, wherein the velocities of initial and final heavy quarkonia v 1.2 are not 
equal to each other, when their product tends to zero, in contrast to the hard limit v± = v<i- 
The spin symmetry relations between the form factors are conserved after the taking into 
account the Coulomb a s /v corrections, which can be written down in covariant form. We 
investigate numerical estimates in the sum rules schemes of spectral density moments and 
Borel transform and show an important role of Coulomb corrections. Next, we perform the 
calculation of gluon condensate contribution to the three-point sum rules of both full QCD 
and NRQCD for the case of three massive quarks, for the first time. 

The paper is organized as follows. The QCD sum rules of three-point correlators are 
considered in Section 2, where the spectral densities are calculated in the bare quark- loop 
approximation and with account of the Coulomb corrections, and the gluon condensate term 
in the Borel transform scheme is presented. Numerical estimates of semileptonic decay modes 
are also given here. Section 3 is devoted to the NRQCD sum rules for recoil close to zero. The 
spin symmetry relations are derived and the gluon condensate is taken into account in the 
scheme of moments. The results are summarized in Section 4. Appendices A and B contain 
technical details in evaluation of decay widths for the massive leptons and gluon condensate 
in full QCD, respectively. 



2 Three-point QCD sum rules 



In this paper we will use the approach of three-point QCD sum rules @, |T(J in the study of 
form factors and decay rates for the transitions B+ — > ip(ri c )l + i>i, where / denotes one of the 
leptons e, \i or r. This procedure is similar to that of two-point sum rules and the information 
from the latter on the coupling of mesons to their currents is required in order to extract the 
values for the form factors. Thus, in our work we will use the meson couplings, defined by the 
following equations: 

(0|OT 5 g 2 |P(p)> = /p f ^ , (l) 

mi + m 2 
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and 

(0\qi^Q2\V{p,e)) =ie fl M v f v , (2) 

where P and V represent the scalar and vector mesons with desired flavour quantum numbers, 
respectively, and mi,m 2 are the quark masses. Now we would like to describe the method 
used. 



2.1 Description of the method 

As we have already said, for the calculation of hadronic matrix elements relevant to the 
semileptonic -B c -decays into the pseudoscalar and vector mesons in the framework of QCD, 
we explore the QCD sum rule method. The hadronic matrix elements for the transition 
.B+ — > ip(r] c )l + i>i can be written down as follows: 

(VMIV^BM) = f + (Pi+P2), + f-q f ,, (3) 
^(J/ii>(p2)\Vn\B c fa)) = iF v e^e w ( Vl +p 2 ) a q^ (4) 

l -(J/^p 2 )\A,\B c ( Pl )) = F A e; + F A (e*-p 1 )(p l +p 2 ) fM + F A (e*-p 1 )q^ (5) 

where g M = (p\ —p-zjp and & = e^(p 2 ) is the polarization vector of J/^-meson. and are 
the flavour changing vector and axial electroweak currents. The form factors f±, Fy, F A and 
F± are functions of q 2 only. It should be noted, that by virtue of transversality of the lepton 
current Z M = ^7^(1 + 75)^7 in the limit vn\ — > 0, the probabilities of semileptonic decays into 
e + z/ e and /i + z/^ are independent of /_ and F A . Thus, in calculation of these particular decay 
modes of B c - meson these form factors can be consistently neglected 0, 0, |T3[. However, since 



the calculation of both the semileptonic decay modes, including e, /i or r, and some hadronic 
decays, stands among the goals of this paper, we will present the results for the complete set 
of form factors given in Eqs. (3)- (5). 

Following the standard procedure for the evaluation of form factors in the framework of 
QCD sum rules, we consider the three-point functions: 

Mvi,V2,q 2 ) = i 2 j dxdye 1 ^*-^ ■ 

(0|T{g 1 (x)75 92 (x),^(0),6(y)7 5 c( 2 /)}|0), (6) 
n^(pi,p 2 ,? 2 ) = i 2 j ' dxdye'to-*-™) ■ 

(0|T{g 1 (x) 7 ,g 2 (a;),J / r' A (0),%)75c( 2 /)}|0), (7) 

where qi(x)j 5 q 2 (x) and qi(x)^ u q 2 (x) are interpolating currents for states with the quantum 
numbers of r\ c and J/ip, correspondingly. J^< A are the currents and of relevance to the 
various cases. 

The Lorentz structures in the correlators can be written down as: 

it, = n+G^+^v + rL^, (8) 

n£, = i^V^ua/3P2Pl, (9) 

u A = n A g„ u + u A p 2 ^ PlyU + n&> 1)MPl> „ + u A p 2 ^p 2tU + u A Pl ^ P2)U . (10) 
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The form factors f±, fy, and F± will be determined, respectively, from the amplitudes 
n±, U v , and = i(II^ ± n^). In (8)-(10) the scalar amplitudes 11* are the functions of 
kinematical invariants, i.e. ITj = q 2 ). 

To calculate the QCD expression for the three-point correlators we employ the Operator 
Product Expansion (OPE) for the T-ordered product of currents in (6)-(7). The vacuum 
correlations of heavy quarks are related to their contribution to the gluon operators. For 
example, for the {QQ) and (QGQ) condensates the heavy quark expansion gives 

(QQ) = — — (G 2 ) L_^(G 3 ) + ... 

YlmQ 7i 360mg Ti 

(QGQ) = ^logK)^(G 2 )--^^(G 3 ) + ... 

Then, in the lowest order for the energy dimension of operators the only nonperturbative 
correction comes from the gluon condensate: 

Il l ( P l,pl,Q 2 )=^ ert (plplq 2 ) + ^\pipl,q 2 )(-G 2 ). (11) 

The leading QCD term is a triangle quark loop diagram, for which we can write down the 
double dispersion representation at q 2 < 0: 

nr rt (p?,#L<? 2 ) = / -A ^V/ 2 '^ \^ ds i ds 2 + subtractions, (12) 

(2tt) 2 J (si -Pi){s 2 -P2) 

where Q 2 = —q 2 > 0. The integration region in (p!2|) is determined by the condition 

_ 1 < 2gis 2 + (gi + s 2 - q 2 ){m 2 b - m 2 c - sj) < 1 
A 1/2 (si, s 2 , q 2 )\ l/2 {m 2 c) si, ml) 

and 

A(xi, x 2 , x 3 ) = (x 1 + x 2 - x 3 ) 2 - Ax x x 2 . 

The calculation of spectral densities pf ert (si, s 2 , Q 2 ) and gluon condensate contribution to 
flTT|) will be considered in underlying sections. Now let us proceed with the physical part 
of three-point sum rules. The connection to hadrons in the framework of QCD sum rules is 
obtained by matching the resulting QCD expressions of current correlators with the spectral 
representation, derived from a double dispersion relation at q 2 < 0. 

H-i(p 2 nPi,q 2 ) = -jttt^ I T 1 ^S'/ 2 ' ~ \ ds i ds 2 + subtractions. (14) 

(Z7T) J [Si — Pi){S 2 — p 2 ) 

Assuming that the dispersion relation (|T4"D is well convergent, the physical spectral functions 
are generally saturated by the lowest lying hadronic states plus a continuum starting at some 
effective thresholds sf 1 and s\ . 

pfv s ( Sl ,s 2 ,Q 2 ) = P r(s 1 ,s 2 ,Q 2 )+ (15) 

e( Sl -sl h )-e(s 2 -sf)- P r nt (s 1 ,s 2 ,Q 2 ), 
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where 



P r(suS 2 ,Q 2 ) = (0|c 7M (75)c|J/^(^))(J/^(r/ c )|F l (g 2 )| J B c )( J B c |6 7 5c|0)) • 

(2tt) 2 (50i - Ml)8(s 2 - Mj) + higher state contributions, (16) 

where Mi )2 denote the masses of quarkonia in the initial and final states. The continuum of 
higher states is modelled by the perturbative absorptive part of IT, i.e. by pi. Then, the 
expressions for the form factors Fi can be derived by equating the representations for the 
three-point functions IT in ([TTJ) and (P~3|), which means the formulation of sum rules. 

2.2 Calculating the spectral densities 

In this section we present the analytical expressions to one loop approximation for the pertur- 
bative spectral functions. We have recalculated their values, already available in the literature 
[f[T |. Among new results there are the expressions for p_,p^ and p±, where p± are spectral 
functions, which come from the double dispersion representation of Tl± = |(n^ ±Hf). These 
spectral densities are not required for the purposes of this paper, but they will be useful for 
calculation of form factors for the transition of I? c -meson into a scalar meson 2 . The procedure 
of evaluating the spectral functions involves the standard use of Cutkosky rules fTB] . There 
is, however, one subtle point in using these rules. At Q 2 < there is no problem in applying 
the Cutkosky rules in order to determine pi(si, s 2 , Q 2 ) and the limits of integration over si, s 2 . 
At Q 2 > 0, which is the physical region, non-Landau-type singularities appear |17], [Tj|, what 
makes the determination of spectral functions to be quite complicated. In our case we restrict 
the region of integration in si and s 2 by s* and s\ , so that at moderate values of Q 2 the 
non-Landau singularities do not contribute to the values of spectral functions. For spectral 
densities Pi{s\, s 2 , Q 2 ) we have the following expressions: 

3 k 

p+(s 1 ,s 2 ,Q 2 ) = ^^{-(Ai + A 2 ) - A;[m 3 (m 3 - mi) +m 3 (m 3 - m 2 )] - 

[2(s 2 A 1 + s 1 A 2 )-u(A 1 + A 2 )] (17) 
■[ml - - + m 1 m 2 - m 2 m 3 - mim 3 ]}, 
3 k 

p-(s u s 2 ,Q 2 ) = -^^{-(A 1 - A 2 ) - k[m 3 (m 3 - mi) - m 3 (m 3 - m 2 )] + 

[2(s 2 A 1 - Sl A 2 )+ M (A 1 - A 2 )] (18) 

11 

■[ml - - + m x m 2 - m 2 m 3 - mim 3 ]}, 
3 

Pv(si, s 2 , Q 2 ) = p^{(2s x A 2 - uAi)(m 3 - m 2 ) 

+ (2s 2 Ai -«A 2 )(m 3 -mi) +m 3 k}, (19) 
p£( Sl ,s 2 ,Q 2 ) = ■^{(m 1 -m 3 )[ml + -(s 1 A 2 2 + s 2 Al-uA 1 A 2 )} 

-m 2 (ml - -y) -mi(m 2 - -y) (20) 



2 The meson at the P-wave level, for which (0|g l i7 jU Q , 2|-P(p)) = ifpP/j,, where P(p) denotes the scalar P-wave 
meson under consideration, and mi ^ m^- 
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+m 3 [ml - -(Ai + A 2 - u) + m 1 m 2 ]}, 

p^(si,s 2 ,Q 2 ) = — ^{m 1 [2s 2 A 1 - wAi + 4A X A 2 + 2A 2 ] 

mim 3 [4s 2 - 2u] + m 2 [2s!A 2 - itAJ - m 3 [2(3s 2 Ai + SiA 2 ) 

-w(3A 2 + Ai) + k + 4AxA 2 + 2A 2 + m 2 (4s 2 - 2u)} (21) 
6 

+-(mi - m 3 )[4sis 2 AiA 2 - «(2s 2 AiA 2 + siA 2 + s 2 A 2 ) 

+2s 2 ( Sl A 2 + s 2 A 2 )]}, 
3 

p A (si,s 2 ,<5 2 ) = —^j^{kum 3 (2m 1 m 3 - 2m\ + u) + I2(m x - m^sjAj + 

fcA 2 [(mi + m 3 )u — 2si(m 2 — m 3 )] + 2A2(A; + 3usi)(mi — m 3 ) 
+Ai[A;M(m 2 - m 3 ) + 2A 2 (A; - 3w 2 )(mi - m 3 )] + (22) 
2s 2 (mi - m 3 )[2£;m 2 - kA x + 3wA 2 - 6wAiA 2 ] - 
2sis 2 (km 3 - 3A2(mi - m 3 ))], 

P ; A ( Sl)S2 ,Q 2 ) = -^{-2(m 1 -m 3 )[(A;-3 MS2 )A 2 + 6 S 2 A 2 ] + 

ku{m\ — m 3 )(2m 2 + A 2 ) + ku 2 m 3 + A 1 [/cM(2m 1 — m 2 — 3m 3 ) 
-2(mi -m 3 )(fcs 2 - A;A 2 + 3m 2 A 2 )] - (23) 
2si[(mi - m 3 ){2kml - QuA x A 2 - SwA^) + 
2s 2 (A;m 3 + 3miA 2 — 3m 3 A 2 ) + A;A 2 (2m 1 — m 2 — 3m 3 )]}, 

p'_ A (si,s 2 ,Q 2 ) = ^{2(m 1 -m 3 )P + 3 MS2 )A 2 + 6s 2 A 2 ] + 

ku{rrii — m 3 )(2m 2 + A 2 ) + ku 2 m 3 + Ai[ku(—2m 1 — m 2 + m 3 ) 

-2 (mi -m 3 )(fcs 2 - £A 2 + 3u 2 A 2 )] + (24) 

2si[(mi - m 3 )(2£;m 3 - 6wAiA 2 + 3wA 2 ) - 

2s 2 (/cm 3 — 3m!A 2 + 3m 3 A 2 ) + A;A 2 (2m! + m 2 — m 3 )]}. 

Here = (si + s 2 + Q 2 ) 2 — 4sis 2 , w = si + s 2 + Q 2 , A x = s\ — m\ + m 3 and A 2 = s 2 — m 2 , + m 2 . 
mi, m 2 and m 3 are the masses of quark flavours relevant to the various decays, see prescriptions 
in Fig. 1. 

We neglect hard O(-) corrections to the triangle diagrams, as they are not available yet. 
Nevertheless, we expect that their contributions are quite small ~ 10% and so, taking into 
account the accuracy of QCD sum rules, the correction will not change drastically our results. 

In expressions ( |T2"D the integration over Si and s 2 is performed in the near-threshold region, 
where instead of a s , the expansion should be done in the parameters (a s /v 13(23)), with 1*13(23) 
meaning the relative velocities of quarks in (be) and (cc) systems. For the heavy quarkonia, 
where the quark velocities are small, these corrections take an essential role (as it is the case for 
two-point sum rules [0, ^0|). The a s /v corrections, caused by the Coulomb-like interaction 
of quarks, are related with the ladder diagrams, shown in Fig. 2. It is well known, that 
an account of these corrections in two-point sum rules numerically leads to a double-triple 
multiplication of Born value of spectral density |23|, |21| . 

Now, let us comment the effect of these corrections in the case of three-point sum rules 
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Figure 1: The triangle diagram, giving the leading perturbative term in the OPE expansion 
of three-point function. 



Q 




Figure 2: The ladder diagram of the Coulomb-like quark interaction. 



13| . Consider, for example, the three-point function n M (pi,p 2 , q) at q = q^ ax , where g^ aa , 
is the maximum invariant mass of the lepton pair in the decay B c — > rjjui. Introduce the 
notations pi = (rn& + m c + Ei, 0) and p<i = (2m c + E 2 , 0). At si — M\ and s 2 — Aff we have 
E\ <C {m,b + ?ric) and E2 <C 2m c . In this kinematics, the quark velocities are small, and, thus, 
the diagram in Fig. 2 may be considered in the nonrelativistic approximation. We will use 
the Coulomb gauge, in which the ladder diagrams with the Coulomb-like gluon exchange are 
dominant. Then the gluon propagator has the form 



(25) 



In this approximation, the nonrelativistic potential of heavy quark interaction in the momen- 
tum representation is given by 

Arc 



V(k) 
b 



i« s (k 2 )g, a s (k 2 ) = bMk2/A2] 
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with rif being the number of flavours, while the fermionic propagators, corresponding to either 
a particle or antiparticle, have the following forms: 



_z(l + 7 o)/2 

Sf(p) 



-z(l-7o)/2 



-ib -^+iO' 
The notations, concerning Fig. 2, are given by 

&i = Pi+i - Pi, h = q%- q%-x, Pn+i = qo=p. 

Integration over p®,p° and q® by means of residues yields the following expression 

00 n dv- 

U^E^q) = Wd^n - w | p /* -V((p i+ t ~ Pif) ■ 

n=Xi=l (27r) -Ei- zO) 

1 1 



(26) 



£ n (9 w iJ qj ' - 

k =xj=x (2vr) 3 (^ - E 2 - zO) \ / 

^((Pn+l - Pn) 2 ) = V(( P - Pn ) 2 ), V(( qi - q ) 2 ) = t/(( qi - p) 2 ), 

where iV c denotes the number of colors, /ii and // 2 are the reduced masses of the (be) and (cc) 
systems, correspondingly. This three-point function may be expressed in terms of the Green's 
functions for the relative motion of heavy quarks in the (be) and (cc) systems in the Coulomb 
field, G®(x,y): 



cgW) = gfn/ - w iJ pt F m 

n=l \fc=l J (27r) 3 (^ - ^ - ZO) 



n-1 

x II ^((P* - P*+i)V PlX_ip " y - (27) 

fe=i 



Comparing the expressions (26) and (|27|) , we find 

n^^,^,^) = 2iV c ^o/Gg i ) (x = 0,p)G| ) (p,y = 0)^. 

= 2N c g^J ' 0^(0, z)G ( §l(z,0)d 3 z. (28) 
For the Green's function we use the representation 

Ei ,y) " £ L5 +1 + 7 (Sj ■ (29) 
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Provided x = 0, only the terms with I = are retained in the sum. Then for the spectral 
density one has 



p fl (E 1 ,E 2 ,q, 



2 > 
max ) 



2N c g^ c 1 (O)v^(O) / *^(p)^ 2 (p 



,c 



,c 



dp 



(30) 



where tyf are the Coulomb wave functions for the (be) or (cc) systems. An analogous expression 
can also be derived in the Born approximation: 



(%(E x ,E 2 ,i 



2 ) 

max I 



2N c g li0 ¥ 1 (0)¥ 2 (0) I ^(p)^ 2 (p 



dp 

(2^)3' 



(31) 



Here stands for the function of free quark motion. Since the continuous spectrum Coulomb 
functions have the same normalization as the free states, we obtain the approximation 



«f(0)tff(0) 



p /1 (E 1} E 2 ,q max ) « p AExi E 2 , q mas , , 

W((0)W£(0j 



Pu( E l, E 2,qmax) C ' 



{ 47ra s 


1 — exp ( 




-r 


47ra s 


\ 3^13 


3f 13 / 




3^23 



1 — exp 



Ana. 



• !3 ) 



(32) 



(33) 



where t>i3, t>23 are relative velocities in the (be) and (cc) systems, respectively. For them we 
have the following expressions: 



^13 



^23 



4mi7TT.3 



Pi — (mi — ^3) 2 



1 



4m 2 m 3 



p\ ~ (m 2 - m 3 ) 



(34) 
(35) 



Eq. (|32]) is exact for the identical quarkonia in the initial and final states of transition under 
consideration. However, if the reduced masses are different, then the overlapping of Coulomb 
functions can deviate from unity, which breaks the exact validity of (|32"D . From a pessimistic 
viewpoint, this relation can serve as the estimate of upper bound on the form factor at zero 
recoil. In reality, this boundary is practically saturated, which means that in sum rules at low 
momenta inside the quarkonia, i.e. in the region of physical resonances, the most essential 
effect comes from the normalization factor C, determined by the Coulomb function at the 
origin. The latter renormalizes the coupling constant in the quark-mesonic vertex from the 
bare value to the "dressed" one. After that, the motion of heavy quarks in the triangle loop 
is very close to that of free quarks. 

In accordance with (||P for the Lorentz decomposition of Pfi(pi,p 2 ,q) we have 



P»(Pi,P2,q 2 ) = (Pi+P2)»P+(q 2 ) +q^P-(q 2 



(36) 



As we have seen, the nonrelativistic expression of p^Ei, E 2 , q^ax) * s proportional to the vector 
(5Vo)> which allows us to isolate the evident combination of form factors f±. The relations 
between the form factors appearing in NRQCD at the recoil momentum close to zero will be 
considered below. Here we stress only that we have 



P+(q 



2 ) 

max) 



p+(q, 



2 ) 

max J 



(37) 
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where the factor C has been specified in (|33|). 

In the case of B c — > J/tplvi transition, one can easily obtain an analogous result for p$(q 2 ) 
(note that the form factor Fq ~ p$ gives the dominant contribution to the width of this decay 
|[L4| ). In the nonrelativistic approximation we have 

PtiQmaJ = Po' B (Qmax) C - ( 38 ) 

To conclude this section, note that the derivation of formulas (^) and (|38|) is purely formal, 
since the spectral densities are not specified at g^ aa . (one can easily show pf(q 2 ) to be singular 
in this point). Therefore, the resultant relations are valid only for q 2 approaching q 2 ^^. 
Unfortunately, this derivation does not give the q 2 dependence of the factor C. But we 
suppose that C does not crucially effect the pole behaviour of the form factors. Therefore, 
the resultant widths of transitions can be treated as the saturated upper bounds in the QCD 
sum rules. 



2.3 Gluon condensate contribution 

In this subsection we will discuss the calculation of Borel transformed Wilson coefficient of the 
gluon condensate operator for the three-point sum rules with arbitrary masses. The technique 
used is the same as in [ITS] with some modifications to simplify the resulting expression. As was 



noted in [jig] , this method does not allow for the subtraction of continuum contributions, which, 
however, only a little change our results as the total contribution of the gluon condensate to 
the three point sum rule is small by itself (< 10%), and, thus, its continuum portion is small, 
too. The form of the obtained expression does not permit us to use the same argument as in 
[8| to argue on the absence of those contributions at all. Their argument was based on an 



expectation, that the typical continuum contribution can show up as incomplete T functions 



in the resulting expression and the absence of them in the final answer leads authors of [ITS]] to 
the conclusion, that such contributions are actually absent in the processes, they considered. 

The gluon condensate contribution to three-point sum rules is given by diagrams, depicted 
in Fig. ||. For calculations we have used the Fock-Schwinger fixed point gauge J22], |20| : 



= 0, (39) 

where A®, a = {1, 2, ...,8} is the gluon field. 

In the evaluation of diagrams in Fig. ^ we encounter integrals of the typef] 

r rl^k h h h 

' C> ~ J (2tt) 4 [k 2 - m 2 } a [( Pl + k) 2 - m 2 ] b [(p 2 + k) 2 - m|] c ' 1 ' 

Continuing to Euclidean space-time and employing the Schwinger representation for propaga- 
tors, 

_ _ = / daa a ~ l e- a{p +m \ (41) 

[p 2 + m 2 } a r(a) Jo v ; 



3 Since the diagrams under consideration do not have UV divergencies, there is no need for a dimensional 
regularization. 
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Figure 3: The gluon condensate contribution to three-point QCD sum rules. The directions 
of pi, ki, k 2 momenta are incoming, and that of p 2 is outgoing. 
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we find the following expression for the scalar integral (n = 0): 

( 1^o+6+c." poo poo poo 

I (a,b,c) = \ ' / / / dadpd^- 1 ^- 1 ^- 1 

r(a)r(6)r(c) Jo Jo Jo 

d k r -a(k 2 +ml)-f3(s 1 +k 2 +2p 1 -k+mi)-7(s 2 +k 2 +2p2-k+mi) 



(2vr) 4 

This representation proves to be very convenient for applying the Borel transformation with 

B p 2(M 2 )e- ap2 = 5(1 - aM 2 ). (43) 

Then, we have 



( — 1 \a+b+c- 



r(a)r(6)r(c)16vr 2 

Z7 (o + 6 + c-4,l-c-6), (44) 
1)a+6+c+1 . , ^ 

r(a)r(6)r(c)i67r 2 ^m 2 c m 2 c 

Z7 (o + 6 + c-5,l-c-6), (45) 
W,cj r(G)r(6) r( c )i 67 r 2 1,M 2 C + M 2 J 

+ l|) ( M ^ 4_a " C ^ 4 " a_6C/ °( a + 6 + c - 6, 1 - c - 6) 

U (a + b + c-6,2- c-b), 

where M 2 C and M 2 C are the Borel parameters in the s\ and S2 channels, respectively. Here we 
have introduced the Uo(a,b) function, which is given by the following expression: 

5-i 



poo 

U Q (a, b) = / dy(y + M 2 C + M 2 C ) a y b exp 
Jo 



Bo - B lV 



(47) 



where 



Bo = jJj^iMUml+m^+MUml + ml)), (4* 



Si 



m 2 



Then, one can express the results of calculation for any diagram in Fig. ^| through Io(a, b, c), 
I^(a,b,c), Ifj, v {a,b,c) and their derivatives over the Borel parameters, using the partial frac- 
tioning of the integrand expression together with the following relation: 

^(^(A^bi] m ^ n w#.(«.M) = 

[Ml]m[MlT d(M 2 bc r d(M 2 c ) n ^ m ^"W( ffl . & ' c )- ( 49 ) 
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The obtained expression can be further written down in terms of three quantities Jo(l, 1, 1), 
1^(1,1,1), 7^(1,1,1) and their derivatives over the Borel parameters and quark masses by 
means of 

T / 1 \ l T / \ / \ 

" lw -^ (o " c) = r(a)r(6)r( c ) rf(mi)- 1 d{m\f^ d(^r^ W./U i, 1. i)- (50) 



However, contrary to the case, discussed in |18|], in such calculations the values of parame- 
ters a, h will arise, for which the Uo(a, 6)-function has no analytical expression (it is connected 
to nonzero m 3 mass in our case). The analytical approximations for Uo(a,h) at these values 
of parameters lead to very cumbersome expressions. The search for the most compact form 
of the final answer leads to the conclusion, that the best decision in this case is to express the 
result in terms of Uo(a, 6)-function at different values of its parameters. For this purpose we 
have used the following transformation properties of Uo(a, b): 



dU (a,b) [ml m\Mf 



^ u ^ b)+ Si u ° (a - b+ih (5i) 

a(7„(a-l,i>)-|-^-^S£](7„( ,i.-l) + 



dU a (a,b) m\ m\M t 

-mr = aU « {a - l ' b) -{m-^L 

2 



^ u ° M+ iik^ b+1) - (52) 

In Appendix B we have presented an analytical expression, obtained in this way, for the Wilson 
coefficient of gluon condensate operator, contributing to the Tlx = Il + + II_ amplitude. One 
can see that, even in this form the obtained results are very cumbersome. So, we have realized 
the gluon condensate corrections as C++ codes, where the functions U (a,b) are evaluated 
numerically. Analytical approximations, which can be made for the Uo(a,b) functions are 
discussed in Appendix B. 

In Fig. |4] we have shown the effect of gluon condensate on the /i(0) form factor in the 
Borel transformed three-point sum rules. 

We can draw the conclusion that the calculation of gluon condensate term in full QCD sum 
rules allows one to enlarge the stability region in the parameter space for the form factors, 
which indicates the reliability of sum rules technique. 



2.4 Numerical results on the form factors 

First, we evaluate the form factors in the scheme of spectral density moments. This scheme is 
not strongly sensitive to the values of threshold energies, determining the region of resonance 
contribution. In the calculations we put 

k th (bc) = 1.5 GeV, 

k th (cc) = 1.2 GeV, (53) 

m b = 4.6 GeV, 

m c = 1.4 GeV, 
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/i(0) 
0.17,— 



0.16 



0.15 

* 

0.14 
0.13 



5 10 15 20 25 30 35 M C 2 C , GeV 2 

Figure 4: The gluon condensate contribution to the /i(0) form factor in the Borel-transformed 
sum rules at fixed M 2 C = 70 GeV 2 . The dashed line represent the bare quark-loop results, 
and the solid curve is the form factor including the gluon condensate term at (-G 2 ) = 10~ 2 
GeV 4 . 

where kth is the momentum of quark motion in the rest frame of quarkonium. The chosen val- 
ues of threshold momenta correspond to the minimal energy of heavy meson pairs in specified 
channels. 

The typical behaviour of form factors in the moment scheme of QCD sum rules is presented 
in Fig. ||. 

The evaluation of Coulomb corrections strongly depends on the appropriate set of ot s for 
the quarkonia under consideration. The corresponding scale of gluon virtuality is determined 
by quite a low value close to the average momentum transfer in the system. So, the expected a s 
is about 0.5. To decrease the uncertainty we consider the contribution of Coulomb rescattering 
in the two-point sum rules giving the leptonic constants of heavy quarkonia. These sum rules 
are quite sensitive to the value of strong coupling constant, as the perturbative contribution 
depends on it linearly. The observed value for the charmonium, « 410 MeV, can be obtained 
in this technique at af ul (cc) = 0.6. The value /b c = 385 MeV, as it is predicted in QCD sum 
rules |]23|| , gives a c ° ul {bc) = 0.45. We present the results of the Coulomb enhancement for the 
form factors in Table [I]. 

The result after the introduction of the Coulomb correction is shown in Fig. [| Such large 
corrections to the form factors should not lead to a confusion, as they are resulted from the 
fare account of the Coulomb corrections both for bare quark loop diagram and meson coupling 
constant. 

In the scheme of the Borel transformation we find a strong dependence on the thresholds of 
continuum contribution. We think that this dependence reflects the influence of contributions 
coming from the excited states. So, the choice of k t h values in the same region as in the scheme 
of spectral density moments results in the form factors, which are approximately 50% greater 
than the predictions in the moments scheme, where the higher excitations numerically are not 
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Figure 5: The QCD sum rule results in the moment scheme for the -^(0) form factor in the 
bare approximation. 




Figure 6: The QCD sum rule results in the moment scheme for the F^(0) form factor with 
account of Coulomb corrections. 
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essential. In this case we can explore the ideology of finite energy sum rules [24|, wherein the 
choice of interval for the quark-hadron duality, expressed by means of sum rules, allows one 
to isolate the contribution of basic states only. So, if we put 



hhibc) = 1.2 GeV, 

k th (cc) = 0.9 GeV, (54) 

then the region of the lowest bound states is taken into account in both channels of initial and 
final states, and the Borel transform scheme leads to the results, which are very close to those 
of moment scheme. The dependence of calculated values on the Borel parameters is presented 
in Figs. [7] and |], in the bare and Coulomb approximations, respectively. 



approx. 


/+ 


f- 


F v , GeV- 1 


F*, GeV 


F£, GeV' 1 


, GeV" 1 


bare 


0.10 


-0.057 


0.016 


0.90 


-0.011 


0.018 


coul. 


0.66 


-0.36 


0.11 


5.9 


-0.074 


0.12 



Table 1: The form factors of decay modes into the heavy quarkonia at q 2 = in the bare 
quark-loop approximation taking into account for the Coulomb correction. 




Figure 7: The Borel transformed sum rule results for the F^(0) form factor in the bare 
approximation of QCD. 

As for the dependence of form factors on q 2 , the consideration of bare quark loop term 
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shows that, say, for F^(q 2 ) it can be approximated by the pole function: 



*3V 



jfr(0) 

1 - 



(55) 



M 2 , 

pole 



with M po i e fa 4.5 GeV. The latter is in a good agreement with the value given in flzfl . However 



we believe that the pole mass can change after the inclusion of a s corrections^ From the naive 
meson dominance model we expect that M po i e ps 6.3 — 6.5 GeV. 




Figure 8: The Borel transformed sum rule results for the F^(0) form factor with account of 
Coulomb corrections. 

We have calculated the total widths of semileptonic decays in the region of M po i e = 4.5 — 6.5 
GeV, which result in the 30% variation of predictions for the modes with the massless leptons 
and more sizable dependence for the modes with the r lepton (see Table [|). 

To compare with other estimates we calculate the width of & — ^ ce + v transition as the 
sum of decays into the pseudoscalar and vector states and findQ BR(_B+ — > cce + u) & 3.4 ± 0.6 
%, which is in a good agreement with the value obtained in potential models and in OPE 
calculations |§, where the following estimate was obtained BR(.B+ — > cce + u) ps 3.8 %. 

In the presented results we have supposed the quark mixing matrix element |Vb c | = 0.040. 



As for the hadronic decays, in the approach of factorization [25| we assume that the width 
of transition — > J/ip(r] c )+ light hadrons can be calculated with the same form factors 
after the introduction of QCD corrections, which can be easily written down as the factor 



4 In HQET, the slope of Isgur-Wise function acquires a valuable correction due to the a s -term. 

5 For normalization of the calculated branching ratios we used the total B c width obtained in the framework 



of the OPE approach [8] . 
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mode 


T, 10- 15 GeV 


BR, % 


i] c e + v e 


11 ± 1 


0.9 ±0.1 


T] C T + U T 


3.3 ±0.9 


0.27 ±0.07 


ipe + v e 


28 ±5 


2.5 ±0.5 




7±2 


0.60 ±0.15 



Table 2: The width of decay modes into the heavy quarkonia and leptonic pair and the 
branching fractions, calculated at tb c = 0.55 ps. 

H = N c a\. The factor a± represent the hard a s corrections to the four fermion weak interaction. 
Numerically we put a\ = 1.2, which yields H « 4.3. So, we find 

BR[B+ -> J/ijj + light hadrons] = 11 ± 2%, (56) 
BR[B+ ^r] c + light hadrons] = 4.0 ± 0.5%. (57) 

Neglecting the decays of b — > ccs, which are suppressed by both the small phase space and 
the negative Pauli interference of decay product with the charmed quark in the initial state 
|| , we evaluate the branching fraction of beauty decays in the total width of B c as 

BR[B+ -> cc + X] = 23 ± 5%, 

which is in agreement with the estimates in other approaches || [|, where this value is equal 
to 25 %. 

3 Three-point NRQCD sum rules 

The formulation of sum rules in NRQCD follows the same lines as in QCD, the only difference 
is the lagrangian, describing strong interactions of heavy quarks. 

3.1 Symmetry of form factors in NRQCD and one-loop approxi- 
mation 

At the recoil momentum close to zero, the heavy quarks in both the initial and final states have 
small relative velocities, so that the dynamics of heavy quarks is essentially nonrelativistic. 
This allows us to use the NRQCD approximation in the study of mesonic form factors. As in 
the case of heavy quark effective theory (HQET), the expansion in the small relative velocities 
to the leading order leads to various relations between the different form factors. Solving these 
relations results in the introduction of an universal form factor (an analogue of the Isgur-Wise 
function) at q 2 -> q 2 max . 

In this subsection we consider the limit 

< + (58) 
W = Vi ■ v 2 — > 1, 
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where t>f 2 — P\,il \jp\,% are the four-velocities of heavy quarkonia in the initial and final 
states. The study of region (|58|) is reasonable enough, because in the rest frame of 5 c -meson 

(Pi — 0))? the four-velocities differ only by a small value \p 2 \ (p% = (yPl,^), whereas 

their scalar product w deviates from unity only due to a term, proportional to the square of 

I t>2 I - w — \ 1 + ~ 1 + ij^J-. Thus, in the linear approximation at Iml — » 0, relations (|58|) 

are valid and take place. 

Here we would like to note, that ([58]) generalizes the investigation of [14], where the case 
of V\ = t>2 was considered. This condition severely restricts the relations of spin symmetry 
for the form factors and, as a consequence, it provides a single connection between the form 
factors. 

As can be seen in Fig. 1, since the antiquark line with the mass m 3 is common to the 
heavy quarkonia, the four-velocity of antiquark can be written down as a linear combination 
of four- velocities v\ and v<i- 

v% = a< + K- (59) 

In the leading order of NRQCD for the kinematical invariants, determining the spin structure 
of quark propagators in the limit w — > 1, we have the following expressions: 



Pi 


-> (mi + m 3 ) 2 = 


•M 2 , 


pI - 


(m 2 + m 3 ) 2 = 


Ml 


A 1 - 


-> 2m 3 M u 




A 2 - 


-> 2m 3 M 2 , 




u - 


2MiM 2 -w, 






-> 4A4 2 A-l 2 (w 2 - 


-I)- 



(60) 



In this kinematics it is an easy task to show, that in (|5Tf) a = b = — ^, i.e. 



T#=-i(«i+t*r. (6i) 

Applying the momentum conservation in the vertices on Fig. 1 we derive the following formulae 
for the four- velocities of quarks with the masses mi and m 2 : 

~< = <+^>i-^ 2 y\ (62) 

v% = v$ + p-(v 2 - Vi y, (63) 
2m 2 

and in the limit w —>■ 1, we have v\ = v\ = 1, as it should be. 

After these definitions have been done, it is straightforward to write down the transition 
form factor for the current J M = Q\Y I1 Q2 with the spin structure = {7^, 757^} 

(H Qi qMHq 2 q 3 ) = tr[T ^(1+^)^(1 +v v 3lu ) ■ 

r 2 ^(l + ^2 A 7A)] •/i(m 1) m 2) m 3 ) ) (64) 
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where r\ determines the spin state in the heavy quarkonium Q1Q3 (in our case it is pseu- 
doscalar, so that r\ = 75), T 2 determines the spin wave function of quarkonium in the final 
state: T 2 = {75,^7^} for the pseudoscalar and vector states, respectively (H = P,V). The 
quantity h is an universal function at w — > 1, independent of the quarkonium spin state. So, 
for the form factors, discussed in our paper, we have 

(P Ql Q 3 \QirQs\PQ 2 Q 3 ) = (cf -V^ + C^V^-h, (65) 

(P Ql Q 3 \QirQs\V Q2 Q 3 ) = ic v -e^e v v la v 2p -h, (66) 
(P Qi q 3 \Qi1^Q3\V Q2 q 3 ) = (c e -ef t + c 1 -v^e-v 1 ) + c 2 -v^e-v 1 ))-h, (67) 

where 



C e = -2, 

m 3 (3m! + m 3 ) 



Cl = 



4mim2 



c 2 = - — (4mim 2 + m x m z + 2m 2 m 3 + ml), 

Amim 2 

cy = — - — - (2m 1 m 2 + mim 3 + m 2 m 3 ), (68) 

2mim 2 
p m 3 m 3 

Ci = 1 + 



= 1 



2mi 2m 2 
m 3 m 3 



2mi 2m 2 

Then for the form factors in NRQCD we have the following symmetry relations: 

f + (c(-M 2 -c^M 1 )-f-(c(-M 2 + c^-M 1 ) = 0, 

Fq • cy — c e ■ FyM.\M. 2 = 0, (69) 

F A (c 1 + c 2 )-c € Mi(F£(M 1 + M 2 ) + F A (M 1 -M 2 )) = 0, 

F A cf + c e -M 1 (f + + f-) = 0. 

Thus, we can claim, that in the approximation of NRQCD, the form factors of weak currents 
responsible for the transitions between two heavy quarkonium states are given in terms of the 
single form factor, say, Fq. The exception is observed for the form factors F+ and F A , since 
the definite value is only taken by their linear combination F+(Mi + M 2 ) + F A (M\ — M 2 ). 
This fact has a simple physical explanation. Indeed, the polarization of vector quarkonium e M 
has two components: the longitudinal term and transverse one e T (i.e. (e*r • t>i) = 0). e£ 
can be decomposed in terms of v± and v 2 : 



where 



ev = 



e» = ael + 04, a 2 + 1 = 1, (70) 



;-2s 2 Pi + up 2 ) -> - I (-< + wv$), 



f S 2 k ' ' y/w 2 - 1 



a = -^^(cl-Pi) -> — , = = (e-^i)- (71) 

'ur — 1 
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From ([n]) one can see that the decomposition of polarization vector e into the longitudinal 
and transverse parts in NRQCD is singular in the limit w — > 1: 

e" = -^t("< + ^)(e ■v 1 )+(3. e£. (72) 
w z — 1 

It means, that the introduction to the form factor Fq of an additional term AF^ = (w 2 — 1) -cm, 
which vanishes at w — > 1 and, thus, is not under control in NRQCD, leads to a finite correction 
for the form factors Fjj 4 and F^. This correction is cancelled in the special linear combination 
of form factors, presented in (|69|). 

In the case of v\ = v 2 we reproduce the single relation between form factors F§ and f±, 
as it was obtained early in fllj . 



Thus, we have obtained the generalized relations due to the spin symmetry of NRQCD 
lagrangian for the case v i ^ v 2 in the limit, where the invariant mass of lepton pair takes its 
maximum value, i.e at the recoil momentum close to zero. 

In the one-loop approximation for the three-point NRQCD sum rules, i.e. in the calculation 
of bare quark loop, the symmetry relations ( |55D take place already for the double spectral 
densities p^ R in the limit \p 2 \ — > 0. We have checked, that the spectral densities of the full 
QCD in the NRQCD limit w — > 1 satisfy the symmetry relations (|69D. 

It is easily seen that in this approximation, 

a,nr 6mim 2 m 3 
\p 2 mi + m 3 ) 



When integrating over the resonance region, we must take into account that 

< 1, (74) 



\p 2 \y/2u 1 m 13 



and so we see, that in the limit \p 2 \ — > the integration region tends to a single point. Here 
Pi = (mi + m 3 + ^i>0), p 2 = (m 2 + m 3 + u 2 ,p 2 ) and m^- = J^^. is the reduced mass of 
system (QiQ 



2 2 

After the substitutions of variables uj\ = tt — and x = (u) 2 — tt — ) • r^fr, in the limit 

1 2mi3 v z 2m23 ' |P2|«' 

\p 2 \ — > for the correlator Uq' nr we have the following expression: 

jjA,NR 1 f duJid£j 2 A,NR _ 

" (27r) 2 J (ujx - wi)(£u 2 — uj 2 ) 

3 /■*«. k 2 dk 



(75) 



where fcth denotes the resonance region boundary. In the method of moments in NRQCD sum 
rules we put u>i = —(mi + m 3 ) + q\ and uj 2 = — (m 2 + m 3 ) + q 2 , so that in the limit qi )2 — > 
we have 



i i jn+m 



n! m! dql'dq™ 



o _ - L - L o 

3 k 2 dk 



TT 2 



(^i + ^ +1 (^ + i) m+1 



(76) 
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In the hadronic part of NRQCD sum rules in the limit \f>2\ ^ we model the resonance 
contribution by the following presentation: 

° ~j-(mi + m 3 )Ml Ml 3 ™ f £ l^MfJ [m^ ) ' [U) 
Saturating the p\ and p 2 channels by ground states of mesons under consideration we have 

Uo' NR [n,m\{ mi +m 3 ) njn 

.pQiQ:i j\/f2 -CQ2Q3 i\/f2 
JlS 1V1 1,1SJ1S 1V1 2,1S 

The value of Fq 1S _^ iS at fixed n = 4 is presented in Fig. |9] at 



0.1S-1S- , Ql Q 3A/r2 fQ,Q 3 irt M l,lS M 2,lS- {'») 



k th = 1.3 GeV, 
m 6 = 4.6 GeV, 
m c = 1.4 GeV. 

Further, in ( |77| ) we can take into account the dominant subleading term, which is the con- 
tribution by the transition of 25 — ► 25. In this case one could expect that the form factor 
is not suppressed in comparison with the contribution by the 15 — > 25 transition, since in 
the potential picture the latter decay has to be neglected because the overlapping between 
the wave functions at zero recoil is close to zero for the states with the different quantum 
numbers []. So, we can easily modify the relation (|T8D due to the second transition and justify 
the value of F^ 2 s^2S ^° reach the stability of F^g^g at low values of moment numbers. We 
find Fq 2S ^ 2S / Fq 1S _^ is w 3.7 and present the behaviour of the form factor F^ at zero recoil 
in Fig. '[Hi 



3.2 Contribution of the gluon condensate 

Following a general formalism for the calculation of gluon condensate contribution in the Fock- 
Schwinger gauge |22|. ^(J we have considered diagrams, depicted in Fig. || In our calculations, 
we have used the NRQCD approximation and analysed the limit, where the invariant mass of 
the lepton pair takes its maximum value. Here we would like to note that the spin structure, 
in the leading order of relative velocity of heavy quarks, does not change, in comparison 
with the bare loop result for the nonrelativistic quarks. Thus, we can conclude, that in this 
approximation, the symmetries of NRQCD lagrangian lead to a universal Wilson coefficient 
for the gluon condensate operator. As a consequence, relations (^) remain valid. 

Below we perform calculations for the form factor Fq(jp\,p\,Q 2 ) in the limit q 2 — > q^ax- 
The contribution of gluon condensate to the corresponding correlator is given by the following 
expression: 

AFf = -^[3^-^, (79) 

7T 4(5 

6 The corresponding estimates were performed in p6| , where the 15 — > 25 transition is suppressed with 
respect to 15 — > 15 as 1/25. 
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Figure 9: The NRQCD sum rule results for the F^(0) form factor. The solid line represents 
the bare quark loop contribution, the short dashed line is the result obtained by taking into 
account the gluon condensate term R only, and the long dashed line is the form factor 
including the full expression for the gluon condensate. 
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Figure 10: The NRQCD sum rule results for the F A (0) form factor. The contribution of 
25' — > 2S transition has been taken into account. The notations are the same as in Fig. [|. 
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where 



Ro 
R2 



k 2 dkdkr> 



ir 2 i J P 3 (k, fc )Pi(fc, fc , Wi)P 2 (fc, k , W2. 



-R, 



a- 



k A dkdko 



n 2 i J P 3 (k,ko)P 1 (k,k ,u 1 )P 2 (k,ko,U2, 
In these expressions the inverse propagators are 

A- 2 



-Rk 



10) 



P 1 {k,k ,LO 1 ) 



U\ + k - 



2mi 
k 2 



P 2 (k,k ,uJ 2 ) = u 2 + k -- , 

2m 2 



P 3 (Mo) 



2m, 



The functions P 9 and P& are symmetric under the permutation of indices 1 and 2 and they 
have the following forms: 



R„ 



Rk 



1 



nil 



1 1 



1 1 

P 2 m 2 



1 1 

P^ + P 

3 1 



1 



P, 2 ' 



3 1 3 1 3 1 

mf Pf + m|P^ + mf PfPg + ^|P 2 3 P 3 

11 11 11 

+ 



+ 



mim 3 P 2 P 3 2 m 2 m 3 P 2 P 2 m\ P 3 4 



52) 



1 



1 



minis P1P3 m 2 m 3 P 2 Pi mim 2 PiP 2 P 2 mim 2 P 2 P 2 



Let us note, that in the calculations of diagrams in Fig. ||] we have used the following vertex 
for the interaction of heavy quark with the gluon 



L v int = -gsKv^K, 

where = ■ ^ and its Fourier-transform in Fock-Schwinger gauge has the form 



After two differentiations of the nonrelativistic propagators 

1 1 



13) 



^4) 



P(p,m Q ) pv + ^ 



15) 



two types of contributions to the gluon condensate correction appear. The first is equal to 

3 



6) 
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and it leads to the term with Rq. The second expression 



v^asG^GlpiekP 



( a sG 2 ) 



k 2 - (v ■ kf 
12 



■<a.O |A| 



12 



17) 



leads to the term with R 2 , which is of the same order in the relative velocity of heavy quarks 
as Ro, but it is suppressed numerically in the region of moderate numbers for the momenta of 
spectral densities, where the non-perturbative contributions (condensates) of higher dimension 
operators are not essential. 

It is easy to show, that the contributions of Ro and R 2 can be obtained by the differentiation 
of two basic integrals: 



E 2 



where k 





r°° k 2 dkdk 


2m 13 2m 23 


TT 2 i 


Jo P 1 P 2 P 3 


ku + k 23 ' 




r°° k 4 dkdk 


2m 13 2m 23 


TT 2 i 


Jo P1P2P3 


ki3 + k 23 



(k 2 13 + k 13 k 23 + k 2 23 ), 
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"23 



R 
R2 



Rg ■ E , 

Rk ■ E 2 , 



where operators R g and Rk can be obtained from R g and Rk in ( |S2]) after the substitutions: 



1 

Pf 
1 

Pirn 
2 

1 

Pi 



-l) n d n 
n\ dcui ' 



-iy 
i\ 



d 



+ 







duo\ duj- 



As a result we have: 



Ro 



R? 



{ 



d 3 



+ 



1 d 3 



+ 



Q3 



+ 



d 3 



6m! duo 3 6m 2 duj 2 2m\ \du 3 duj 2 duo 2 



2m 2 
3 



{ 



' d 3 
d 4 



+ 



9 3 



A\m{ dto\ 



3\m 2 , 



' d 4 
duo 



+ 



1 (&_ 

Am 2 m 3 \duj 2 



du 2 duji 

_3 d*_ 

4!m| duj 2 
d 4 N 
du 2 diJi / 
d 4 

+ 2 



+ 



2m 13 2m 23 
ki3 + k 23 ' 
3 / d 4 



+ 



+ 



3!mf \du 4 du) 3 duj 2 



+ 



+ 



1 (&_ 

Amiin 3 \dujf 
d 4 \ 



+ 



4!m? 



+ 4 



du 2 duji duj\duj 2/ 
d 4 „ <9 4 



+ 6- 



9c; 4 duo 3 duo 2 duj 2 duo 2 



+ 4 



+ 2 



d 1 



+ 



d 1 



duj 3 duj 2 duj[d<jj 2 



dijj\dijj\ 



+ 



d 4 ' 
9a; 4 



+ 



*9) 



(90) 
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a 4 a 4 ' 

dujfdu)2 dio\dLo\ dio\dio\ i 

QA QA Qi « 

+ 3— t— r + 



+ ^ ) ~ 1 f ^ I 

2mim 2 V9o;f9a;2 ' du^du^ duidu^ I 4mim 2 \du 2 duj2 



2m 13 2m 2 3 
mi3 + m 23 



(kf 3 + k 13 k 23 + k? 
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Equations ([75|), fl5D| ) and flSO]) represent the most compact analytical expression for the con- 
tribution of gluon condensate to the form factor Fq, whereas performing the differentiations 
leads to very cumbersome expressions. 

In the moment scheme of sum rules we suppose 



L0 X 
L0 2 



= -(mi+m 3 ) + gi, 
= -(m 2 + m 3 ) + q 2 , 



and expand AF{f 2 in a series over {gi, q 2 } at the point {0,0}, which allows us to determine 
AF^ [n,m] = fq™dq m • Further analysis has been performed numerically with the 

help of MATHEMATICA. The evaluation of the AF^ 2 [n, m] dependence in a broad range of 
[n, m] takes too much calculation time, so we restrict ourselves by showing the results in Figs. 



and 10 for the fixed n = 4 with the following set of parameters 



m b = 4.6 GeV, 

m c = 1.4 GeV, 

-G\ v ) = 1.7 ■ 10- 2 GeV 4 . 

7T M 



, a 



by taking into account (Fig. [|) and without accounting (Fig. |10]) the contribution of transition 
2S — > 2S, where F Q (2S — > 2S)/F (1S — > IS) ~ 3.7. As can be seen in these figures, the gluon 
contribution, while varying the moment number in the region of bound cc-states, plays an 
important role, because it allows us to extend the stability region for the form factor F^ up 
to three times (from m < 5 till m < 15), and, thus, the reliability of sum rule predictions. 
Let us also note, that in the scheme of saturation for the hadronic part of sum rules by the 
ground states in both variables sj and s 2 , the account for the gluon condensate leads to the 
20% reduction for the value of form factor F^. 

The analysis of the dependence on the moment number m in the region of bound states 
be at fixed m shows that the contribution of gluon condensate in this case does not affect the 
character of this dependence. This may be explained by the fact that the Coulomb corrections 
to the Wilson coefficient of gluon operator G 2 ^ play an essential roleQ. The summation of a s /v- 
terms for be may give a sizeable effect, contrary to the situation with the cc system, where the 
relative velocity of heavy quarks is not too small. 

7 At present, there is an analytical calculation of initial six moments for the Wilson coefficient of gluon 
condensate in the second order of a s pTj. In this region, the influence of gluon condensate to the sum rule 
results is negligibly small for the heavy quarkonia, containing the 6-quark, which does not allow one to draw 
definite conclusions on the role of such a s -corrections. 
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To conclude, we have calculated the contribution of gluon condensate to the three-point 
sum rules for heavy quarkonia in the leading order of relative velocity of heavy quarks and 
in the first order of a s . Due to the symmetry of NRQCD in the limit, when the invariant 
mass of the lepton pair takes its maximum value, i.e. at the recoil momentum close to zero, 
the Wilson coefficient for the form factor Fq is universal in the sense that it determines the 
contributions of gluon condensate to other form factors, in accordance with relations (|69|). 

4 Conclusions 

We have calculated the semileptonic decays of B c meson in the framework of sum rules in 
QCD and NRQCD. We have extended the previous evaluations in QCD to the case of massive 
leptons: the complete set of double spectral densities in the bare quark-loop approximation 
have been presented. The analysis in the sum rule schemes of density moments and Borel 
transform has been performed and consistent results have been obtained. 

We have taken into account the gluon condensate contribution for the form factors of 
semipeltonic transitions between the heavy quarkonia in the Borel transform sum rules of 
QCD, wherein the analitycal expressions for the case of three nonzero masses of quarks have 
been presented. 

We have considered the soft limit on the form factors in NRQCD at the recoil momentum 
close to zero, which has allowed one to derive the generalized relations due to the spin sym- 
metry of effective lagrangian. The relations have shown a good agreement with the numerical 
estimates in full QCD, which means the corrections in both the relative velocities of heavy 
quarks inside the quarkonia and the inverse heavy quark masses to be small within the accu- 
racy of the method. Next, we have presented the analytical results on the gluon condensate 
term in the NRQCD sum rules within the moments scheme. 

In both the QCD and NRQCD sum rules, the account for the gluon condensate has allowed 
one to enforce the reliability of predictions, since the region of physical stability for the form 
factors evaluated has significantly expanded in comparison with the leading order calculations 
of bare quark-loop contribution. 

Next, we have investigated the role played by the Coulomb a s /v corrections for the semilep- 
tonic transitions between the heavy quarkonia. We have shown that as in the case of two-point 
sum rules, the three-point spectral densities are enhanced due to the Coulomb renormalization 
of quark-meson vertices. 

The complete analysis shows that the numerical estimates of various branching fractions: 



agree with the results obtained in the framework of potential models and Operator Product 
Expansion in NRQCD. More detailed results are presented in tables. 

Thus, we draw the conclusion that at present the theoretical predictions on the semileptonic 
decays of B c meson give consistent and reliable results. 
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5 Appendix A 

In this Appendix we present the derivation of exclusive semileptonic widths for the 5 c -meson 
decays into the J/ip, rj c mesons with account of lepton masses. 

The exclusive semileptonic width Tsl for the decay B c — > J/ip(rj c )Wi, where I — e, \i or r, 



can be written down in the form 2c 



v "=ih^Tnr 1**1 (91) 



where d 4 q = 2n\q\dq 2 dq , dr\ = \f>i\dQ,i / (IQn 2 \fcp) is the leptonic pair phase space, dVLi is the 
solid angle of charged lepton /, \p{\ = y/q I §i/2 is its momentum in the dilepton center of mass 



(]9T ) are given by 



system and = Jl — 2A + + XL, with A± = (mf ±m^)/g 2 . The tensors L af3 and W a p in Eq. 



U* = \ E {lO a vi){viOPl) = 2|pfpJ + p?p° - g a \ V i ■ P VI ) + ie a Wpi 7 p„ lS }, (92) 



spins 



and 



Wop = J d ^\ Pl -p 2 -q)W a(L 



d A p 2 

6{E 2 )8{Ml - 2Mi -q + q 2 - M 2 )W a p\ P2=Pl . q , (93) 



where 



W af} = (f + (t)( Pl +p 2 ) a + f-(t)Qa)(f + (t)(Pi+P2)p + (94) 
for the pseudoscalar particle in the final state, and 

W a/3 = -{iF^g ae + iF+p le (px + p 2 ) a + iF^p le q a - F v e aei j{p x + p 2 )Y) • 
(iF A + iF^ Plk ( Pl + p 2 ) p + iF^p lk q p + F v ep k mn(pi + p 2 ) m g n ) • 
( E ^) (95) 

polarizations 

for the vector meson in the final state with the polarization e M , where 

E ^ = j^-9 ek ( 96 ) 

polarizations 2 

M 2 is the mass of final state meson (J/ip or rj c ). 
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The integral over the leptonic phase space in Eq.(^T|) is given by 



dT l L°^ = ^-^(L aP ), (97) 



with 



[L a ?) = ^ / dJW* = |{(1 + AOCgV - <? Q V) + ^A 2 <fV}, (98) 

where Ai = A + — 2A 2 and A 2 = A + — A?_. Introducing the dimensionless kinematical variable 
i = q 2 /ml and integrating over g the semileptonic width takes the following form: 

1 G 2 AVh l 2 m? /■*«»»<« a ~ 

where 

1 



2 

In Eq. (|99|) the limits of integration are given by t m i n = and t max = — ?(Mi — M 2 ) 2 . 
Calculation of (L al3 )W a p yields the following expressions 

(L aP )W aP = 1 -(3t 2 f 2 _(t)X 2 mi + 6tUt)f + (t)(M 2 -Ml)X 2 m 2 b + 

/ 2 (t)(t 2 (2Ai - 3A 2 + 2)m\ - 2t(M 2 + M 2 )(2Ai - 3A 2 + 2)m 2 b + 
2(M 2 - M 2 ) 2 (Ax + 1))) 

for the pseudoscalar meson in the final state, and 

{L af3 )W aP = -^-(2(t 2 (A 1 + lK-2t((l + A 1 )M 2 -5(l + A 1 )M 2 + 9M 2 A 2 )m 2 + 

(M 2 - M 2 ) 2 (l + Ai))^) 2 - 2(t 2 m\ - 2t(M 2 + M 2 )m 2 + 
{Ml - M 2 ) 2 )(F^(2(1 + Ai)im 2 - 3£A 2 m 2 - 2(1 + Ai)M 2 + 
2(1 + Ai)M|) - 3tF A m 2 A 2 )F,f + (t 2 m\ - 2t{M 2 + M 2 )m 2 b + 
(M 2 - M 2 ) 2 )((2t 2 m 4 b + 2t 2 X 1 mt - 3t 2 A 2 m^ - 4(1 + Ai)tM 2 m 2 - 
4(1 + \ x )tMl + 6tM 2 A 2 m 2 + 6tM 2 X 2 m 2 b + 2(1 + \\)M* + 
2(1 + Ai)M 2 4 - 4(1 + Ai)M 2 M 2 )(F^) 2 + 
6tF^F^X 2 m 2 b (M 2 - M 2 ) + tm 2 b (3t(F* ) 2 X 2 m 2 b + 
16(1 + Ai)(F y ) 2 M 2 2 - 24A 2 M 2 (F y ) 2 ))) 

for the vector meson in the final state. 



6 Appendix B 

In this Appendix we illustrate the kind of expressions, which arise for the gluon condensate 
contribution to the form factors F l (t) in the framework of Borel transformed three point sum 
rules, in the case of = f+(t) + f-(t) form factor. 
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Following the algorithm, described in section 2.3, for H[ we have the following expression 

nf 2> = c{- 1 '- 2) c/ (-i,-2) + ci- 1 '- 1) c/ (-i,-i)+ E c[°' i) u (o,i) + 

i=-4 

E Cf'^o(M) + E Cf^o(2,i) + E C?'^o(3,i), (100) 



where 



i=— 5 i=— 6 i=— 6 





- 2 ) = 


m 2 c + m; 






12M 2 C 




-i) = 


1 






12M 2 ' 


cf- 


-4) = 


1 ( 

48M fe 6 c M c V 



(M™(6m 2 + ml) + 8m 2 2 M* c M* c + (-5m 2 + 5m 2 + 



-m 2 + 2mim 2 - m^ + 4m 2 m 3 - 16m 2 + M 2 C + Q 2 )M 2 c M t 



2 

cc ' 



(101) 
(102) 



m\)Mt c Mt c + (-7m 2 + 3m 2 + 2 mi (m 2 - m 3 ))M 4 M C 6 C - 4m 2 M 2 M C 8 C (103) 
-2m 2 M c 1 c °), 

Cl(0 '~ 3) = 4* J> M* ((m ' ~ 2mim3 ~ 2m 3) M ce + (6m 2 + 4M 2 c )M fe 8 c - 
(2m 2 + 2m 2 , + 6m 2 m 3 + 2m 2 + 2mi(— 5m 2 + m 3 ) — 11M 2 C 
-2Q 2 )M^ C M 2 C + (m\ - m 2 + 2mi(m 2 - 5m 3 ) - 2m 2 m 3 + 2m 2 (104) 
+5M 2 C + Q 2 )M b 2 c M c 6 c + (-7m 2 + 4m!m 2 - 2m 2 - lOm^ - 
8m 2 m 3 + 2m 2 + 12M 2 C + 2Q 2 )M 4 M 4 C ), 

C ''~ 2) = 48^ M 4 ((3m ' ~ 2mim3 ~ 2m * )M ~ + (m2 ~ 5m ' + 9 ^)^ + 
(—4m 2 — 7m 2 , + 2mi(6m 2 — 5m 3 ) — 8m 2 m 3 — 13m 2 + 16M 2 C + 
4Q 2 )M fe 4 M 2 C - (5m 2 - 4m!(m 2 - 4m 3 ) + 2(m 2 + 2m 2 m 3 + m 2 - (105) 
Mi~Q 2 )MlMi), 

Cf _1) = 4gM 6 M 4 (-2^M 4 + (-6m 2 + 5M 2 C )M 4 + (106) 



( o,o) m 2 (M 2 c + M 2 c ) 
Cl " 48MPS ' (107) 

C ' {1 '" 5) = ~ 48M* M6 (m2M ' c ~ m i M ce)(( 5 ^2 + m§)M£ - m 2 M b 6 c M c 2 c ~ (108) 
(4m 2 + 3m 2 )M 4 M 4 + 2m x (m 2 - m 3 )M 2 c M c 6 c + 2m 2 M 4 ), 

°^ 4) = 48M* M6 (m ' (3m ' + 2mim3 + 6m 3) M c C ° + 

(-2m 4 - 7m 2 M 2 c + m^M^jM^ + (8m 4 + 4m 3 m 3 + 3m 2 m 2 + 

2m 2 m 3 + m\{Am\ + m\) — 2mi(4m 2 — 3mlm 3 + m 2 m 2 — m| — (109) 

llm 2 M 2 c - 3m 2 M 2 c - 4m 2 Q 2 - mjQ 2 )M* c M 2 cc - (3m 4 + 6m 2 m 2 + 
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L cc 



m\{— 61712 + 4m 3 ) + 2mi(m 2 m 3 + 2m 2 M 2 c — 2m 3 M 2 c ) + 
m^lOm 2 . + 2m 2 m 3 + 2m\ - 8M C 2 C - 3Q 2 ))M 4 M C 6 C + (km x m\m 3 + 
m\(9rn\ + 16M 2 C ) + m 2 2 (2m 2 2 + 4m 2 M - 3 - 4m 2 . - M 2 C - 
2Q 2 ))M fe 6 c M 4 c - mi (6m? + 4m 2 m 3 + 2m 2 (-m 2 + m 3 ) + + 
2m 2 m 3 - 4m 2 + 5M 2 C - Q 2 ))M 2 C M C 8 C ), 

^ 1,_3) = 48M 8 M6 ( 2m 3(^i + 5m 2 m 3 + 2m^)M c 8 e + (m 4 + m 2 (9m 2 - 4M 2 e 

8m 2 m 3 M 2 - 3m 2 M 2 c - 11M 4 )M 8 + (2m 4 + 10m 2 m 2 + m\ + 
7m 2 M 2 c + 18m 2 m 3 M 2 c + 2m 2 M 2 c - 20M 2 C - 4 mi (m 2 - m 3 )(m 2 (110) 
+2M 2 C ) + m 2 (2m 2 + 5M 2 C ) - 2m 2 Q 2 - 5M 2 c Q 2 )M fe 6 c M 2 c + (-7m 2 m 2 
-4m 4 + 4m 2 M 2 c + 10m 2 m 3 M 2 c - 7M 4 + 2m 2 (2m 2 m 3 - m 2 + 5M 2 C ) + 
m 1 (Qm 2 - 6m 3 - 2m 2 m 3 + 4m 2 m 2 + 4m 2 M c 2 + 10m 3 M 2 c ) + 2m 3 Q 2 
-4M 2 C Q 2 )M 4 Ml + (m 4 - 2m\{m 2 - 4m 3 ) + 10m 2 M 2 c + 2m 1 m 3 (m 2 + 
m 2 m 3 + m 2 - M 2 C - Q 2 ) + m 2 (m 2 + 2m 2 m 3 + 7M 2 C - Q 2 )), 

°^ 2) = 48 M6 ( m '(3m 2 - 2m 1 m 3 + 6m 2 )M c 6 c + (14m 2 m 3 M 2 c + (16m 2 - 

be cc 

UMDMl + m 2 (4m 2 + M 2 c ))M fe 6 c + (-5m 2 m 2 - 2m 2 m^ - 2m 4 + 

4m 2 M 2 c + 8m 2 m 3 M 2 c + 32m 2 M 2 c - 8M 4 + m 2 (-3m 2 + 4M 2 C ) + (111) 

2m 1 (3m 2 m 2 - 2m 3 - 2m 2 M 2 c + 6m 3 M 2 c ) + 3m 2 Q 2 - 4M 2 C Q 2 )M 4 M 2 C 

— (2m' 1 m 3 + m 1 (— 4m 2 m 3 + 4m 3 + M cc ) + m 3 (m 2 + 2m 2 m 3 + 4m 3 

-7M 2 C - g 2 ) + 2m 1 m 3 (m 2 - m 2 m 3 + 3m 2 - 2M 2 C - Q 2 ))M 2 C M C 6 C ), 

= ~ 48M« M6 ( 2mim 'M c 4 c + (4m 4 - 17m 2 M 2 c + M 4 )M 4 + 

2m 2 (m 2 - 2mim 2 + m 2 + 2mim 3 + 4m 2 - 6M 2 C - Q 2 )), (112) 
r (2,-6) (M 2 c -M 2 c )(m 2 M 4 -m 2 M 4 ) 3 

Cl " 96M« ' (U3) 

Cf^ = ~ 96M 8 M e (^M 4 - m 2 M 4 )((m 2 + 2m 2 )M b 6 c + 7m 2 M 2 c M 4 + 

m 2 (m 2 - 2m 1 m 2 + m\ + 2mim 3 + 2m 2 m 3 + 3M 2 C - Q 2 )M 4 C - (114) 
m 1 (m 1 — 2mim 2 + m 2 + 2mim 3 + 2m 2 m 3 + 7M CC — Q )M CC ), 

Cl(2 '~ 4) = 96M^ ^ m \ m l M l! + ml(6m 2 M 2 c - 6m 3 M 2 c + r^m 2 - 

2m 1 M 2 c - 4m 3 M 2 jM 2 c M 8 - (12m 3 2 m 3 M 2 cc + ±m 2 m\M 2 cc - 12m 2 M 4 
-4m 2 M 4 + m 4 (3m 2 + 4M 2 ))M 1 C ° + 2m 1 (-3m 2 ; M 2 + 4m 1 m 2 m 3 M 2 
-Sr^M 4 + m 2 (3m 3 M 2 c - m x m\ + 3 mi M 2 cc + 2m 3 M 2 c ))M b 6 c M 4 + 
2m 2 (-4m 2 (m 1 - 2m 3 )M 2 c + m\{ra\ + 3M 2 C ) + (3m 2 + 6)71^3 + (115) 
2M 2 C - 3Q 2 )M 2 C )M 4 Ml + m 2 (m 2 (m 2 - 4M 2 C ) + 4m 2 (m 1 - 4m 3 )M 2 
+2(-2m 2 - 4 mi m 3 + 3M 2 C + 2Q 2 )M 2 cc )Ml M 2 C , 
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Cf ' 3) = 48M 8 M 6 K(9m 3 M 2 c - r^m 2 + m 2 (-9M 2 c - 4 mi m 3 + 4m 2 ) + 

4miM c 2 )M b 2 c M c 4 c - (2m^m 3 + 9m^ + m 4 + 10m 2 m 3 M 2 c - (116) 
10M 4 )M b 6 c + (m 4 m 2 + 2m^m^ - 3m 2 m 2 M 2 c - 3m 2 M 4 - 14m 2 m 3 M 4 
+9M C 6 C + m 2 (m 2 m 2 + 2m 2 M 2 c - 3M 4 C ) + 2mi(-m^m 2 + m\m\ + 
m 2 m 3 M c 2 - 2m 2 m 2 M 2 c + m 2 M c 4 - 3m 3 M 4 ) - m 2 m 2 Q 2 - 2m 2 M 2 c g 2 + 
3M 2 c g 2 )M b 4 + m^mfm 2 + 3(m 2 - m 3 )m 2 M fe 2 c + m 2 (3m 2 M b 2 c - 3m 3 M fe 2 
-2m 2 m 2 + 2m 3 3 ~ 2m 3 M cc) + ^i^ 3 (m 2 + 2m 2 m 3 - 3M 2 C - g 2 ))M c 4 ), 

^ = 96M io M 8 (^ m l m t M cc + m^MlMi (6m 2 M 2 - 6m 3 M 2 + r^m 2 - 

2mi M 2 + 4m 3 M 2 C ) + M fe 6 c (3m 2 m 4 + 8m 2 m 3 (m 2 - M 2 C ) M 2 C + (117) 
4M 4 (-11m 2 + 2M 2 C )) + M 4 m 3 M 2 c (-8m 2 (m 1 - 2m 3 )m 3 M 2 + 
m 2 (m^ + 6m 3 M 2 c ) - 2M 2 c (-3m 2 m 3 - Qm^l + 2m l M 2 cc + 
20m 3 M 2 +3m 3 g 2 ))), 

Cf -1) = 96M 8 M e m 2 (M 2 c (4m 2 m 3 + 15m 2 - 16M 2 C ) + m 3 (-m 2 m 3 + (118) 

2mi(m 2 m 3 - m 2 + 2M 2 C ) + m 3 (-m 2 - 2m 2 m 3 + 13M 2 C + Q 2 ))), 
r (2,o) rn 4 (M 2 c (m 2 -4M 2 c )+m 2 M 2 c 

96M£°M« ' lU9j 

(3,-6) _ (M 4 m 2 - m 2 M 2 c ) 3 
Cl " 96« ' (120) 
C[ 3 >~ 5) = Q6M s M 8 (M 4 m 2 - mjM^(Mtmi(2m 2 m 3 + M 2 e ) - 

m 2 M 4 (2m 2 m 3 + 7M 2 c )), (121) 
^ 3 '" 4) = ~ 96M^°M8 ( ~ 3m ' m 3^ c 8 c + M b 8 c m 2 (m 2 m 2 - 8m 2 m 3 M cc + 2M 4 ) + 

2M 4 m 2 M 2 c (m 2 m 2 + 6m 2 m 3 M 2 c + 3M 4 )), (122) 
^ 3 " 3) = 48M« M« K^'M 4 (-2m 2 m 3 + M 2 ) + M^ c (—2m\m\ + 2m 2 m 2 M 2 c + 

6m 2 m 3 M 4 -3M C 6 C )), (123) 
(3,-2) _ m 2 (3m 2 m 2 M c 4 + M 4 (m 2 m 2 + 12m 2 m 3 M 2 c - 18M 4 )) 

r (s,-i) m|(2m 2 M - 3 - 9M 2 ) 

Cl " ' (125) 

r (3,0) _ m 3 / 12fi x 

Cl " 96M-MB- (126) 
For functions U (a,b) we have the following expressions 

C/oM) = 1 E +fe 2Q^- 1 exp[-5 ](M 2 + M 2 c ) a+6+1 - 2 K_ n [2^ B^B-,}, 

n=l+b \ Bi J 
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for a > 0. Here K n [z] is the modified Bessel function of the second order. In the case 
of U (— 1, — 2) and U (— 1, — 1) we have failed to obtain exact analytical expressions, so we 
present their analytical approximations: 



U (-l,-2) = 5^{2exp[-^]^ + 2exp[-<Bf 2] ^^^ 



+ exp[-B 1 _/ 2 B 3 1 /2 ]B 3 _ /2 B 1 1 /2 (2 + vB l )v - 



2v 2 B_ lB f exp[^-i]r 0, ^ + + F(0, B^ 2 ) 



-2, J B_ 1 S 1 (exp[^i]r (o, ^ + ^ j + r(o, <5f 2 ) 

- exp[t;5 1 ]r(0, Bi^S^ 2 + vB 1 )) - v 2 B_ x exp[£ V2 £f 2 ] (127) 
—B 2 _ l B 1 {v 2 B 2 Y{Q, B]/ 2 B 3/2 ) + 2uSir(0, B]! 2 Bl 12 ) 
+2r(0, B]! 2 b\ /2 - 2 exp[uBi]r(0, B]! 2 b\ /2 + uSi))}, 

v \ v \ B ± J 
T(0, B l _! 2 Bl /2 ) - expt» J Bir(0, B^Bl 12 + uSi) + 

,(exp^r(0,^ + ^)-r(0,^)) 5l ) (128) 



-y^ wex p[- s -i 2s i 8/a ] + 5-i(^ir(o,B V2 Bf 2 ) + 

T(0, fl^flf/ 2 ) - exp^^O, B^Bf 2 + ^))}, 

where t> = M^ c + M 2 C and r(a, z) is the incomplete gamma function, which is given by the 
integral T(a, z) = / 2 °° t a ' l e^dt. 
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